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into account, thus warranting fulfilment of the optical theorem. The relevance of this analysis to 
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I. INTRODUCTION 



In recent years there has been a rapidly growing interest in the field of particle optics, especially neutron and 
atom optics (for a recent review see and Q respectively, and references quoted therein), due to a spectacular 

improvement of the experimental techniques, connected to the introduction of the single crystal interferometer in the 
first case, and to progress in microfabrication technology and development of intense tunable lasers in the second one. 
Such new achievements provide very important tests verifying the validity of quantum mechanics, especially in that 
it predicts wavelike behaviors even for single microsystems. 

At the same time a new challenge arises, linked to the accuracy required in the description of the interaction 
between the microsystem and the apparatus acting as optical device. The question of the description of the dynamics 
of a microsystem interacting with a system having many degrees of freedom (e.g., matter seen as an optical medium 
characterized by an index of refraction) has been extensively studied and contains some typical quantum mechanical 
features, such as quantum correlations between the two systems, by which a reduced description of the microsystem's 
degrees of freedom can arise only by suitable approximations. This subtle point is particularly important in the 
case of particle optics, where the main interest is devoted to the coherent wavelike behavior of particles, as can be 
justified on the basis of the similarity between a Schrodinger equation with an optical potential and the Helmholtz 
wave equation |,|. The very existence of such an optical description of the interaction is far from trivial and strongly 
depends on the experimental conditions. The attention has been mostly devoted to exploiting the optical analogies, 
while little has been said on the borderline between the optical regime, in which coherent effects are predominant 
and a classical wavelike description plays a major role, and an incoherent regime, where incoherent effects, caused by 
the interaction between the microsystem and the apparatus and showing typical particlelike features, should not be 
neglected. This attitude is exemplified in neutron optics by the use of the "coherent wave" formalism, instead of a 
reduced density matrix description, as usually adopted in quantum optics. 

In this paper we want to address the question of how to consistently describe both regimes applying a recently 
developed approach to the description of irreversible subdynamics in quantum mechanics to the specific case 

of neutron-matter interaction. In this approach the use of an effective T-matrix describing the local interactions 
as practical starting point leads to the introduction of a time scale and in the particular case of particle-matter 
interaction to a dynamical semigroup, whose generator has the typical Lindblad form | p^ . The expressions appearing 
in the generator are linked to particle-particle interactions, like the Fermi pseudopotential, and to properties of the 
macroscopic system, like the dynamic structure function, first introduced by van Hove [ pl| . The first part of the 
generator accounts for the description of the coherent interaction in terms of optical potential and index of refraction 
well-known in neutron optics P,0,13[ . The remaining part is shown to be related to the dynamic structure function 
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or, equivalently, to the density correlation function and leads in a straightforward way to results obtained in the so 
called "rigorous theory of dispersion" Q . 

The paper is organized as follows: in Sec. II we give an account of the formalism; in Sec. Ill it is applied to neutron 
optics; in Sec. IV we consider diffuse scattering, connection to the dynamic structure function and fulfilment of the 
optical theorem; in Sec. V we evaluate possible experimental consequences; in Sec. VI we comment on our results 
indicating potential future developments. 

II. INTRODUCTION OF THE FORMALISM 

In this section we briefly introduce the formal scheme restricted to the description of a microsystem following , 
to which we refer the reader for further details. We indicate by Ti'^^ the Hilbert space in which the microsystem 
is to be described; its energy eigenvalues are Ef, with energy eigenstates uj, spanning Ti,^^\ Both systems will be 
considered confined, e.g., in a box. We shall adopt the second quantization formalism, setting for the Hamiltonian H 
of the system: 

f 

where a/ is the destruction operator for the microsystem, either a Fermi or a Bose particle, in the state w/; is 
the Hamilton operator for the sole macrosystem {[H^.af] = 0). Indicating by Hp the whole Fock space and by H° 
its subspace in which iV = aj^cih^ the number of microsystems, is equal to zero, we will denote with |A) the basis 
of eigenstates of iJm spanning Ti", -ffm|A) — E\\X), N\X) = 0. V represents the interaction potential between the two 
systems. Having in mind to describe situations in which only one particle is observed in each experimental run, or 
equivalently a collection of noninteracting particles in each run, we assume for the statistical operator the following 
expression: 

9f 

where is a statistical operator in the subspace W° , representing the macrosystem, and therefore 

a^e"" = e"'a} = V/, 

while (3 is a statistical operator in the subspace Ti.^ of Tip in which iV = 1. The coefficients ggf build a positive, trace 
one matrix, which can be considered as the representative of a statistical operator g in TL^^\ Being interested in the 
subdynamics of the microsystem we shall exploit the following reduction formula, valid for any operator of the form 

Trup {Ag) = ^ ^/ggg/ = Tr^d, (^Agj . 
f-.g 

Wc wish to determine the equation driving the time evolution of the statistical operator on a time scale t much longer 
than the typical duration of microphysical interactions for the macrosystem, and therefore we shall approximate 
by: 

Trnp (^a}age-^"^g{t)e^"^^ - ggfit) . 

To proceed further we will exploit the cyclicity of the trace operation, shifting the time evolution on the destruction 
and creation operators, thus working in Heisenberg picture. In this way no simplifying assumption is made on the 
structure of p™. We now introduce the following superoperators, that is to say mappings acting on the algebra 
generated by creation and destruction operators: 

h' = ^[H,-], n'o^^[Ho + Hn,,-], V'^^[V,-]. (2.1) 

Making use of this mappings we evaluate (aJiOfc) with the aid of the following integral representation: 



[ggf{t + T) - ggf{t)] 
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+ioo+r; 



dz 
27ri 



a-k, 



Let us stress at this point the relevance of the formaUsm of second quantization. The operator quantities of interest 
can be expressed in terms of products of creation and destruction operators. The study of their time evolution may 
thus be reconducted to evaluation of field operators of the form '^a\ connecting in the Fock space subspaces with n 
and n + 1 particles (and similarly for '^a^ connecting subspaces with n and n — 1 particles). Thus, even recovering 
at the end the usual one particle quantum mechanics, the Fock space structure plays a central role and accounts for 
the similarities between this simple case and the description of macroscopic systems . For the mappings defined 
in (2.1) identities hold that are reminiscent of the usual ones in scattering theory: 



{z-n') '^{z-n'o) [i + v'{z-H') 

In particular we can introduce the superoperator T{z) 



l + iz 



{z-n',)- 



satisfying 



and 



{z-wy 



r(z) = v' + v'(z-H')"'v', 



[z-n',) ' + {z-H'o)-'T{z){z-H'oy 



(2.2) 



(2.3) 



T{z) = v' + v'iz-n'o)-'nz), 



(2.4) 



corresponding to the Lippman-Schwinger equation for the T-matrix. Taking into account the fact that [H, A^] = 
one can see that the restriction to Til. of the operator T[z)ak has the simple general form: 



ihT{z)ak\-Hi^ ^^T'I {ihz) 

h 



(2.5) 



where (z) is an operator in the subspace 7i^. This restriction is the only part of interest to us, since we are 
considering a single microsystem. Our formalism points to this matrix, whose entries are operators on the Hilbert 
space of the macrosystem, as the basic mathematical tool to describe the physics of the microsystem: we will show that 
it yields all relevant quantities and, in our opinion, could be a sou nd s tarting point for phenomenological assumptions. 
TJJ {z) bears a connection to scattering theory, as it is clear from (2.4); it is also related to the thermodynamics of the 
macrosystem being an operator on TiPp. To help clarifying this connection we consider a simple case in which (z) 
can be explicitly calculated. Let the macrosystem be composed of free particles: 



H„ 



V 



E 



where is the creation operator of a particle in an eigenstate Vr, with energy E^, (either a Bose or a Fermi particle) . 
Reca lling that we are describing a single particle and exploiting the superoperators introduced in (pT|), (2.2) and 
(2.3) we can calculate (z) as defined by ( |2.5| ). To do this we bring to normal order the creation and destruction 
operators associated with the macrosystem and restrict ourselves to a one mode dynamics, in which, apart from 
statistical corrections, only one creation and one destruction operator of the type b appear: that is to say we neglect 
three particle collisions. Then one obtains: 



T)(i?fe+*£) = ^fet^fc,^|F(2)+F(2) -VL\f,v)h,. 



where e is a positive quantity and the following relationships hold. 



(fc,^|i/(''|/,^) 

{Ki\VL\f,v) 



{Ef + E^)SkfS^n 
(l±6j&5)(fc,e|V^(''|/,??) 



il±blb^)Vka^; 



(2.6) 



here the superscript (2) denotes operators in the two-particle Hilbert space and statistical corrections for scattering 
in the medium are taken into account in the potential term Vl, implicitly defined by ( 2.6D and by the usual resolvent 
series (+, — sign stand for Bose and Fermi statistics respectively). The connection to the familiar T-matrix is evident. 
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We now come to the master equation describing the irreversible time evolution of the statistical operator on the 
chosen time scale: 



dgkh 
dT 



{Ek - Eh) 6kh-^Yl QkfQfh + QkgQlg + tY1 (^A?)fc/ Sfa^Lx£)lg, 



(2.7) 



/ 



/s 



from which we can read off the structure of the generator of the semigroup driving the time evolution. The quantities 
appearing in (2.7) are defined in the following way: 

Q,/ =Tr„, [T){Ek+ie)Q'-\T)\ 

{\\T){Ek+iemt)) 



Ek+Ex-Ef-E^~ie' 



(2., 



with e a positive constant and ^(r) a complete system of eigenvectors of f?™(T) with eigenvalues 7r^(r)- If we now 
introduce in H^^^ the operators Hq, Q, Laj and g, 



{g\Ho\f) = EfSgf, {9\Q\f)=Qgf, (glU^lf) = (L^^) {dlQlf) = Qgf 



eq. (p^) becomes: 



dgjr) 
dT 



Ho + Heff, gir)] - i {f , gir)} + i ^ U^girp^^ 



where 



H 



Q + Qt 



off 



r = i 



Q-Qt 



Verification of the conservation of the trace of the statistical operator within the adopted approximations leads to the 
following relationship 



and therefore to: 



(2.9) 



dgjr) 

dT 



AC 



(2.10) 



This master-equation is a typical result of the formalism restricted to the case of a single microsystem, for the general 
structure see [|[^] 

Before applying (2.10) to a concrete physical situation it can be useful to gain some further insight into the structure 
of the operators appearing in it. As already said the quantity that the formalism suggests as a natural candidate 
where to put in suitable phenomenological expressions is the operator (z), an operator whose trace over the Fock 
space for the macrosystem calculated with g'^ gives the value of the T-matrix for scattering from state Uf to state Uk 
averaged over the state of the macroscopic system. A quite general phenomenological expression may be obtained in 
the following way. Suppose that T{z) has the form 



T{z) = - [Vizhz), ■] , Viz) - J2 VkXf^i^)alb{afb^, 

kXffi 



with b\ b creation and destruction operators in the Fock space for the macrosystem. We thus have 

ihT{z)ak = '^Vk\ft^{ifiz)b\afb^^ = ^ Ty {ifiz)af, 
A/M / 
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and supposing translation invariance in the interaction kernel: 



(2.11) 



A/1 



Such an Ansatz amounts to introduce an effecti ve po tential which should give in Born approximation the full scattering 
amplitude. As a result the potential term in ( 2.10 ) is linked to the scattering amplitude, as we shall see in the next 
paragraph, while the incoher ent contribution is generally connected to the scattering cross section. To realize this let 
us consider the last term of (^!7|), keeping the proposed Ansatz into account: 



t{Ek +ie,x- y) 



Ex - ie 



Uf(y){\"\i,Hx)i^{x)\\){\\g^^{T)\>!) 



xgf,{r) I dV / dy{\'\^-^{x')i,{x')\\")u;iy') ^^y'), 



(2.12) 



and let us specialize to the case of a diagonal matrix element. Supposing the statistical operator for the microsystem 
is quasi-diagonal and the macrosystem is at equilibrium, so that g'^\\) — p™|A), we exploit the usual representation 
for the delta function, thus obtaining: 

27r 



j:j:^iSiE, + Ex-Ef-Ex, 



f AA' 



d-^x / d'^y ul{y){X\i;' {x)t{Ek + ie,x - y)i;{x)\\')uf{y) 



In this formula one has the typical transition probability between an initial state /, A' and a final state fc. A, averaged 
over all possible initial configurations and summed over all possible final states for the macrosystem, that is to say 
contributions from both coherent and diffuse scattering are included. It might be instructive to show in a different 
way the connection between the last term of (2.7) and the total scattering cross section, referring to a famous paper 
by van Hove pl| ]. Taking for concreteness the Fermi pseudopotential (see next paragraph), whose Fourier transform 
is simply the constant V = ^^j-^, we evaluate the diagonal element of ( ^.12 ) assuming that the Uf are given by plane 
waves (the indexes f,g,h,k becoming momenta), thus obtaining [N{x) = ip'' {x)ip{x)]: 



A, A' 



d^p r d\ 

{2nhf J {2nhf 



d^x / d^y 



{X'\N{x)\\) 



Eu+Ey-^iP+\q) -Ex-ie 



{P+lq\g\P-lq) 



X (A|e-(r)|A) 



2m V ' 2 

{X\N{y)\X) 



Ek + Ey-^{P-^q) -Ex+^s 



and supposing g such that the energies in the denominators may be considered approximately equal, introducing the 
Wigner function for the neutron 



fwix,p) 



one easily has 



n \2TThj 



d^P I dt I d^e-^i^^-'^'-y+rA^-'-'^-r 



X Jd^XU{X,P){N (x-0iv(x + ^,t)), 



(2.13) 



where (...)= Tr-^p (. . . p™), and N{x,t) denotes the operator in the Heisenberg picture. We have thus recovered the 
typical factorized structure appearing in the expression for the scattering cross section of a neutron off a macroscopic 
system: square modulus of the Fourier transform of the interaction potential times the dynamic structure function 



5 



depending on transferred momentum and energy, with the refinement that it is here weighted according to position 
and momentum distribution of the incoming par ticle . For the non-diagonal matrix element one can expect to obtain 
analogous results if the quantities appearing in ( 2.12| ) are sufficiently slowly varying functions of their arguments, so 



that, in the continuous limit, an interpolation formula of the form 



with g(^) a suitably smooth function may be used. The failure of such an approximation and thus the relevance of the 
actual value of the parameter e in the final expression might be traced back to t he breakdown of the approximations 
that have led to the markovian evolution generated by the master-equation ( ^.7| ). 

III. OPTICAL BEHAVIOR 

We now devote our attention to the interaction of neutrons with matter. This field is well suited to test our 
formalism both because of the very refined experiments that have been carried out in neutron interferometry [p]jl^] 
and because of the very well studied description of neutron optics phenomena, as developed for example in the book 
by Sears [Q, that we will take as basic reference. As a first step we want to consider the coherent interaction of 
neutrons with matter and therefore we neglect in ( p.7D the last contribution, linked to incoherent processes. As we 
will see later this term implies indeed a smaller correction in the case of neutron scattering. We are left with: 

= -T {Ek - Eh) Qkh -tY^ QkfQfh + -^Y1 ^i^aQhg: (3-1) 
^ fa 

and we need a suitable expression for the operator 

Qkf = TTH. [T'}iEk+ie)g^{r)]. 

Following Sears we adopt the Fermi pseudopotential to describe the neutron nucleus interaction in impulse approxi- 
mation; let us recall the form of the T-matrix in the context of the elementary theory of dispersion: 

r=^^5„^53(x-R,), (3.2) 

a i—1 

where x is the position operator for the neutron, Ri the position operator for the i-th nucleus of type a, ha the bound 
scattering length, depending on isotope and spin orientation, m the neutron mass, the number of nuclei of type 



a. An operator of the form (3.2), that is to say a sum over one-particle operators, is expressed in second quantization 
by: 

a •' 

where is the field operator, acting in the Fock space of the macrosystcm, corresponding to particles of type 

a. For the sake of simplicity from now on we will consider one kind of particles, thus dropping the subscript a. 
Furthermore we will assume that 6 is a real quantity, since we are not going to deal with absorption phenomena. 
As we shall see in the next section we concentrate on non-hermiticity of the potenti al c onnected with incoherent 
processes and not with net absorption. A phenomenological description as given by ( |3.3| ) falls within the class of 
effective potentials considered in the previous paragraph and corresponds to the following interaction kernel: 

t{z,x-y) = bS^x-y), (3.4) 

m 

leading to 

T'}{Ek+ie) = ^bjd^x J d^y^\x)ul{y)6^x~y)uf{y)^{x). 



Eq. (p.lD thus becomes, in operator form: 



6 



dt 



Hq, g{T) 



i 2Trh 



i 2T:h 



fi m 



b / d^a; (^t(a;)V;(a;))^ ^(t)5^(x - a;), 



(3.5) 



where x is the position operator for the neutron and {A)^ = TiUp. {g ™{T )A). If we consider only pure states and 
assume the macrosystem to be at equihbrium ( {■ ■ = (■••) )j cI- (pj) is equivalent to the following stationary 
Schrodinger equation 





2m 



2TTh^ 



b{^^x)^{x))\cj,{x)^Ec^{x), 



(3.6) 



which, remembering that the average particle density 5'^{x — Ri)) is given in second quantization by {ip'' {x)iIj{x)) , 
is exactly the equation used by Sears to describe all coherent neutron optical phenomena, here recovered in a straight- 
forward, alternative way, though in a very different framework. The term 



2-Kh' 



5(V'+(a;)V'(a;)) 



is called optical potential and assumes different expressions according to the structure of the system. If the medium 
can be considered homogeneous, with density rio, eq. (3.6) describes propagation of matter waves with an index of 
refraction given by 



n=\l- 



2Trh^ 
mE 



brio 



1 - —brio 
Zn 



(3.7) 



as first obtained by Goldberger and Seitz ||T^ in the absence of absorption. This is the formula currently used to 
calculate phase shifts in neutron interferometry experiments [Q: 



(3.8) 



where D is the thickness of the sample. 

In a similar way we can obtain from (3.1) a more general formula for the refractive index introduced for the first 
time by Lax Starting from the general expression (2.11) the potential term in (3.1) becomes 

^QkfQfh{T)^Y^TTH^ d^x / d^yij''{x)ul{y)t{Ek+ie,x-y)uf{y)ij{x)g''\T)gfh{T). 
f f •' •' 

Following Lax we suppose that the system is homogeneous, so that 

Tr„, [^:'{x)^:{x)g^{T)] = n^. 

We have 

^Qkfgfh{T) "^rio^ d^xt{Ek+ie,x) j d^y ul{y)uf{y)gfh{T) 

= no J d^xt{Ek +ie,x)gkh{T), 

where we have exploited the orthogonality between the states {w/}, thus obtaining the matrix element of the T- 
operator for forward scattering, averaged over the possible states of the macrosystem. Keeping the relation between 
T-operator and scattering amplitude into account we come to 

-no f(0,Ek)gkh[T). 

m 



Inserted in the Schrodinger equation this term is equivalent to an index of refraction of the form: 



mEk 



nof{0,Ek 



l + —nof{0,Ek), 
Ztt 



(3.9) 
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simply linked to the forward scattering amplitude. An analogous result holds for electromagnetic waves propagating in 
a material with low density A similar treatment has been proposed |jl^ and adopted (see for example in the 
description of the propagation of atoms through a dilute medium, showing the interest of similar descriptions also for 
atom optics. In the case of thermal neutrons the scattering amplitude is isotropic within a very goo d ap proximation 
and is given in terms of the scattering length by the simple formula f = —b which reduces (3.9) to (3.7). 

So far we have shown how, starting from (2/7) and neglecting the incoherent term, we can recover some important 
results obtained within the framework of multiple scattering theory and used to describe the coherent interaction of 
neutrons with matter. Our formalism puts in evidence the statistical operator of the macrosystem, the T-matrix and 
the scattering amplitude, so that phenomenological inputs are rather direct. Further improvements of the formulas 
obt aine d are allowed by the presence of £'™(t) and depend on its evaluation. The correction factor c that Lax includes 
in (3.9) to obtain the index of refraction 



n - 1 + —nocf{0,Ek), 

2,71 



is connected to fulfilment of the optical theorem, which in our formalism, as we will see in the next section, is related 
to the presence of the incoherent contribution. 



IV. THE INCOHERENT CONTRIBUTION 



We come now to the main statement of this paper, the connection between the contributions other than the 



commutator in (2.10) and the dynamic structure function, toge ther wit h th e relevance of this relationship to the 



optical theorem. As observed by Sears an expression of the form (3/7) or (|3.9| ) for the refractive index doesn't include 
the contribution to the attenuation of the coherent wave in the medium due to diffuse scattering and, hence, violates 
the "optical theorem" of scattering theory To overcome this difhculty he refrains from ad hoc assumptions as 

in p9[| , which amount to introduce a suitable im agin ary contribution to the potential, and considers a rigorous theory 



of dispersion. In this more accurate treatment (3.2) is replaced by 

and /q has the general expression (k is the incident neutron momentum) 

/« = + +0(fc2), 

where the second term had been previously omitted because of its smallness, since typically j^kb < 10^*. Furthermore 
the scattering amplitude is to be multiplied by a constant c which should take local field corrections into account and 
whose value depends only on the temperature, density and chemical composition of the medium. Sears obtains an 
estimate for this constant in terms of the structure function of the macroscopic scatterer in the case of an homogeneous 
medium, applying a multiple wave formalism to solve the scattering problem, and drawing strong analogies to the usual 
descriptions of propagation of electromagnetic waves. In this way he recovers a correspondence between attenuation 
of the coherent wave in the medium and diffuse scattering. In the following we shall set fa = f ^ct and consider 
only real b, in order to concentrate upon diffuse scattering, neglecting absorption. By diffuse scattering we intend all 
scattering that is not coherent in the absolute sense, that is elastic and coherent (for the distinction between absolute 
and relative incoherence see for example P,12|). To compare with these more refined results we have to consider all 



contributions in (2.10). Let us stress from the very beginning some general features of this expression, thanks to 
which it can describe more general physical situations than those arising in an evolution driven by a Schrodinger-like 
equation. The last two terms 

^ELl^?.4+iELAC^L^ (4.1) 

allow for the presence of a non-self-adjoint potential which is nevertheless not linked to real absorption. This is 
the case for the present treatment, in which the imaginary part of the optical potential is to be traced back to the 
existence of diffuse scattering, as opposed to the coherent wavelike behavior. Attenuation of the "coherent wave" 
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is due to the presence of the anticommutator term, responsible for the imaginary potential, balanced by the last 
contribution, typically incoherent in that it leads from a pure state to a mixture. This last term is given by a sum 
over subcollections, formally similar to the expression that we would obtain for the statistical operator after the 
measurement of a given observable (see Q). The subcollections are denoted by the indexes A^, which specify a change 
of the state of the macroscopic system, caused by interaction with the microsystem, thus making this contribution to 
the dynamics incoherent. In fact wc will sec in the case of neutron-matter interaction that the trace of this term gives 
all the contrib utio ns to incoherent scattering, that is to say the total diffusion cross section. The balance between the 
two terms of (4.1) accounts for fulfilment of the optical theorem. 



To see this let us now consider (4.1) in more detail. Starting from ( |2.8| ) and (2.11), introducing a Laplace transform 
for the energy dependence of the effective T-matrix 



t{E,x) 



together with the following expression for the density number operator in terms of creation and destruction operators 
with specified momentum 



N{x) = ijj^{x)^{x) 



V 



.bp-i 



we obtain 



da / d^x" e-i''e-^^"^^-''h{a,x" -x)e^^°^e-^'^-''" 



where V is the volume of the region in which the system is supposed to be confined. Indicating by t{E, k) the Fourier 
transform of the potential with respect to space 

iiE,K) = J d^xt{E,x)e-i'^-'' 

and after some simple manipulations one comes to 

^ 2 

dTe-^UiHo + ie,K)e^^^e^'^^e-^'^-^(\\e-^"^-^b^ 



.bp- 



to be inserted in (4.1). Before doing this let us introduce the useful notation 

e-i"-^Aei"-^ = ^ e-i^\A)A, {A)a =J2\E + A){E + A\A\E){E\, {A)\ = {A')^a. 



We have 



,,P k',P' 



xe R^'^gef 



dr e n'^e'i 2m'^e'i ™ i(Ho + ie, k) 



E 

A. A' 



e-^^^Tr«, 



.bp-.] 

A 



b'* ,&p,^.' 



-A' 



and similarly for the anticommutator part. An important simplification takes place if one can use symmetry under 
time and space translations. Time translation invariance occurs if, at least with reference to the interaction with the 
microsystem, matter can be considered at equilibrium, then: 

Tin, (AAg"S_A') = Saa'T^tu, [Aaq^'B-a) ■ 
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Similarly space translation invariance implies 



Tr-H^ 



-A 



such a symmetry can be implemented at equilibrium in the thermodynamic limit and can be practically assumed for 
a microsystem interacting with a homogeneous portion of a macrosystem. Then one has, performing also the t, r' 
integrals: 



n 1 2 



s 



El 



^ + ^ - A - i£ 



X t(Ho + ?£, n) 



1 



^ + ^ - A + 



-e fi' 



^ + ^ - A + 

m 2m 



iilP + ^ _ A - i£ 



where 



-A 



or equivalently, introducing the x, p dependent amplitude 



in the form: 



-^E( l^'^^^'^o + *^''^'^) 

K.A \ ^ 



(^-^-aV 

\ m Ini J 



-i(Ho + ie, K,x) 



1 



2ee '^'^''tMo + i£,n,x) — 5 p 

^ ' 1^ - ^ - A + i£ 

m zrn 



1 



^ - - A - ze 

m 2m 



-F(Ho + ie, /«,x)e^'^-« /^(k, A). 



(4.2) 



(4.3) 



Introducing this explicit representation in ( 2.10 ) one gets the typical master equation of Brownian motion, that can 
be further simplified in the assumption of small momentum transfer, i.e., expanding the expression with respect to 
K ■ X and K ■ p. Exploiting the fact that £'^^(0, A) contains a (5a,o factor, one can immediately see by inspection that 
the K = contributions cancel each other provided the effective T-matrix is a slow function of energy, 

(fel^l f)t{Ek,0)i\Ef,0) « {k\g\f)^[i\Ek,0)iiEk,0)+PiEf,0)i{Ef,0)] ; 

on the other hand for a homogeneous medium the n = contributions are equal to those obtained writing the 
correlation function as a factorized product 
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Tr 



-A 



provided we assume the condition of "normal density fluctuations", {{N^) — {N)'^)/V'^ -C n^- Instead of restricting 
the sum to the k ^ contributions we can therefore subtract from the correlation function its factorized part. After 
straightforward manipulations, using 

E&kf^^-t = / d'x^^{x)i.ix)ei---, Trn, [(A) ^g"^ {B) _ ^] = J ^e"^^* {BA{t)) 



we come to 



Tr-Hf 



Tr 



(^E^P+t^P-f) ^''"(E^P-t^^+f) 



P-* I Q 

' A 



Tr- 



■Hf 



E^^ 



where 



and finally: 



^''~^'^*J '^^^ J '^^ye^'"^''^^^ {5N{y)5N{x,t)) 
{N{x)) = Tin^ [N{x)q^] , 6N{x) = N{x) - {N{x)) , 



1 I 1 

n 



UELULA«,4 + iELAC^L^ 

K,A \ ^ m ' 2m 



A - Z£ 

1 



it (Ho- 



le.K] 



X t(Ho + ie, — ; ^ 

m 2m 



1 



^ + - A + i£ 



i(Ho + i£, K)e n'^x^ 

1 



^ + A - 

m zm 



-^^^ j d?x j (fye^"- '' {SN{y)6N{x + y,t)) . 



(4.4) 



Thanks to the last term of ( ^.10 ) it is possible to take into account collisions that modify the state of the macroscopic 
system (see 0). The probability per unit time of such coUisions is given by the trace of ^ \ La^pLa^ as seen in 
Sec. II. In the case considered this trace may be written as 



2n rio 



(fk / (fn {k \g\ k) \i{Ek,K - k) 



(4.5) 



X dt (Tx 



j3„^-i(f;:-^)«+i('^-fc)-3; 



Sy-{5N{y)5N{x + y,t)). 



thus recovering again the van Hove structure for the scattering cross section [compare ( 2.13| )], with the difference that 
now the system is considered to be homogeneous, so that only the momentum distribution of the incoming microsystem 
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is of relevance. Let us observe that subtraction of the uncorrelated part of the response function accounts for the fact 
that only diffuse scattering, that is scattering that does not leave the macroscopic system unchanged contributes 
to this term. We now specialize to the case of neutrons, adopting the Fermi pseudopotential given by (3^), so that 
(4.5) becomes 



J J d\ J d^K {k \g\ k) Q - fe] , i 



2m 



2m 



where, denoting by lu and q energy and momentum transfer respectively, 

Sciq,u;) ^^JdtJ d^a.e-'C-'-'J-) J d'y {6N{y)5N{x + y,t)) ^ 



(4.6) 



(4.7) 



If the momentum distribution o f the incoming particle is suitably peaked around po with respect to the momentum 
dependence of 5c we have from ( |4.6| ) 



Tim? 



J d\ j d^K {k \q\ k) Sc Q [po - fe] , ^po - = J Q [po ~ fc] , ^pa 



in particular, in the static limit expression (4^) becomes: 



where 



,2^0 f ^r, a r \ P° 
Hob — / dn„Sc(q) = no—(7d 
m J m 

Sc{q) = d'xe'1-^J d\ {5N{y)5N{x + y)). 



and we have denoted by q the momentum transfer and by the total diffusion cross section per particle. This is the 
result derived by Sears for the attenuation of the coherent beam due to incoherent scattering, which he obtains by an 
evaluation of the local field effects, neglected in the equation giving the optical neutron dynamics (3.6), see |^,||,|l8|. 
In our approach, however, the incoherent contribution is already present in the equation giving the dynamics of the 
microsystem, being connected to the thermodynamic properties of the macrosystem through the response function 
Sc{q,i-u)- This new feature is obtained by means of the more general formalism adopted, leading to a master-equation 
of the Lindblad type for the statistical operator, in which due to the optical theorem a close correlation exists between 
the incoherent contribution and the imaginary part of the optical potential which is not connected to absorption. To 
see this correction to the optical potential let us exploit the simple relation 



A^At, B = Bt, 0=A + i:B 



\Jg- g\J^ = 



A, g 



+ i 



{B,g} 



and write the commutator and anticommutator term of (2.1C) in the form — j- (^Ug — ^0^^ The calculation of is 



essentially given by the anticommutator at the r.h.s. of (4.4) and the commutator in (3.5). In the case of the Fermi 
pseudopotential, using (^) one has: 



U = n„ 



or in the static limit 



An 



^ 2nn' 

U = Uo 



b~i— d\ \k){k\ / dLo^ dn^-sJ-[K-k],- 



2m 



2m 



b-i— f d\ \k)(k\ - I dn„Sc(q) 



(4.8) 



(4.9) 



where q denotes as usual the momentum transfer. Neglecting diffuse scattering we would have U — riob, simply 
a c-number giving the usual refractive index; the remaining part is, in a sense, induced by the optical theorem. To 
compare with the results derived by Sears we have to consider the expression obtained for the static limit ( |4.9| ) applied 
to a plane wave of momentum pq, which gives an idealized description of the preparation of the incoming microsystem, 
thus leading to 
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U = Ho 



(4.10) 

Ch. 4), however 
). The analysis 

that we have put forward rehes on the assumption that the main contribution to the dynamics is given by the 
commutator term in (2.1C), while the terms in (4.1) may, as a first approximation, be neglected. This leads to an 



this expression agrees with the results obtained relying on the idea of local field corrections (see Js 
here ([1.9[) is a direct consequence of the equation driving the dynamics and of the Ansatz (2.11 



optical description, as for the case of neutrons, in which, considering the dimensionless parameter ^^^n^b, the terms 
other than the commutator are of second order. The opposite situation takes place if the interaction is such that the 
main contribution is given by (1.1), while the commutator may be neglected. This happens when dissipative effects 
are predominant, as in the case of Brownian motion mentioned below (4.3), where incoherent interactions through 
collisions involving energy and momentum transfer play the main role, a case we intend to deal with in a future paper. 



V. EXPERIMENTAL IMPLICATIONS 



We now address our attention to potential experimental implications of the above introduced description of neutron- 
matter interaction. Of course possible new features in the dynamics are linked to the presence of the last two terms 
m the r.h.s. of ( ^.lOD , as given by mA) , and such corrections will be generally small, being of second order in ^Uob 

(typically n-ob < 10^^ at thermal neutron energies). In this respect interferometric experiments, in which the 
experimental setup is conceived in order to enhance the coherent behavior, should be particularly relevant: think 
for example of the beautiful experiments realized by the Ranch group in Wien exploiting the perfect crystal neutron 
interferometer [pU^, [l^. 



Consider now eq.( 2.1C ): the map on the r.h.s. is afhne and trace preserving, and therefore clearly predicts neutron 
conservation. Nevertheless the last contribution which offsets the anticommutator term is linked to diffuse scattering: 
one has neutron conservation if also diffuse particles contribute to the experimental observation. This is not so for 
interferometric experiments. In such cases only the wavelike behavior affects the observed dynamics, and thus only 
the commutator part of the evolution map is of relevance: the net result is an imaginary correction to the coherent 
scattering length as in ( |no| ), that is to say a reduction of the neutron flux responsible for the interference pattern. 
This fact is usually taken into account adding an imaginary part proportional to the total scattering cross-section 

|2C|]) according 



to the phase shift calculated as in (3 
to the formula: 

X = X +iX 



thus including both absorption and diffuse scattering (see 



-UobXD 



D 

iriocrt — , 



exp(ix) = exp —iuobXD 



D 

riocrt — 



In the absence of absorption this correction is considered negligible and the relevant incident flux is often evaluated 
simply closing one of the two beam paths. This attitude is however at least in principle incorrect, as it appears taking 
the whole dynamics as given by (2.1C) into account. In fact when one closes the path without the sample also diffuse 
neutrons, which are lost for the interference pattern, having their path "labeled" by scattering with the sample, may 
contribute to the transmitted intensity. The experimental device is no more acting as an interferometer and therefore 
cannot select only those neutrons that have undergone coherent interactions. This additi onal contribution to the 
transmitted neutron flux is given by the trace of the last term of ( 2.10 ), that is to say by (4^). In calculating the 



amplitude of the interference pattern one should therefore rely not simply on the measured transmitted flux, but on 
this quantity minus the additional incoherent contribution given by (4.£ )i thus obtaining a reduction of this amplitude: 
the purely "optical" treatment leads in principle to an overestimate of the visibility of the interference pattern. This 
is normally not the case in real experiments, since the angle of acceptance of diffuse neutrons is very small, as for the 
perfect crystal neutron inte rferometer. Let us give some quantitative estimate of the aforementioned eff ect. 

In order to evaluate (4.6) we have to make a definite choice for the structure function Sc{q, oj), in fact (4.6) is given 
by: 



A 



1, 

h 



-Tr 



„(i) ^ LAfpLt^j 



d\ / d^K {k\q\k)Sc 



k] 



1 




fc2- 




n 


2m 


2m 


) 



where the quantity A takes diffusion at any angle into account, 
isotropic medium, such as a liquid or a gas, one has [H: 



In the static approximation, for a homogeneous and 



5c(g,c^) = S'c(g)5(u;) 



5c(g) = l + noydVe''?'-[g(r)-l], 



(5.1) 
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where g{r) is the pair correlation function. A possible choice for g{r), allowing Sdq) to be evaluated analitically, is 
the following, valid for a dilute hard sphere gas with atomic diameter a: 



9ir) = 



r < a 

1 r > a 



The quantity of interest for us is A in its dependence from the maximal angular acceptance if, determined by the 
experimental apparatus, multiplied by the time the neutron takes to go through the sample. Supposing the momentum 
distribution of the incoming particle sufficiently well peaked around Pq we rewrite A introducing the expression given 
by (^^) and multiplying by the time interval, thus coming to: 



Ai^p) = 2Trnob^D / d6'sin6' 
Jo 



2'Knr,a 



(1 



h 

apo 



3in(if-y2(] 



^v/2(l-cos( 



where cos9 = {pq ■ k)/pfy The primitive of this integral can be straightforwardly evaluated by a change of variables, 
and exploiting the fact that in our model iS'c(O) = 1 — |7ra'^no, we have an explicit representation of diffuse scattering 
at any angle Lp: 



A{ip) = 2Trnob'^D { (1 - cos^) + 3[1 - S'c(O)] 



\aPo 



sin \^^y'2{l-cos(p) 
^V2(l-cos=^ 



considering in particular small ip the expression may be approximated as: 

A{p) ^ nn^b'D f^p'SM + ¥''(2^11 - ^c(O)] (^)' - 3^^0(0)) + 0(^«) 

Let us now consider the experiments performed using the perfect crystal interferometer. The angular acceptance 
is very small, only a few microradians for thermal neutrons pl| ]. Taking for instance a gaseous sample, an order of 
magnitude estimate gives A{ip) ~ 10~^*, that is to say an extremely small quantity, in agreement with the accuracy 
obtained using this interferometer based on Bragg diffraction. An interferometer based on a different physical principle 
could possibly lead to a higher angular acceptance, thus enhancing this effect connected to diffusion. In view of the 
next equation (5.2) a completely different situation arises if one considers systems with abnormally large density 



fluctuations, as would be the case near a first order phase-transition. 

Another point of interest is the linear dependence on 5*0(0) of the leading term in A{ip). The quantity 6*0(0) is 
particularly relevant from the physical point of view, being connected to the isothermal compressibility xt and to the 
fluctuations in the number of particles in the sample 22 1: 



S*c(0) = UokBTxT 



(AiV)2 



N 



(5.2) 



The actual value of S*c(0) cannot be measured experimentally from scattering experiments, and has to be obtained 
by an analytical continuation. The analysis we propose could provide an independent way to measure S*c at q = 0. 
In fact in the static approximation, independently of the particular form of Sc{q), for very small |q|, that is to say for 
very small ip, one has in good approximation: 

A{p) ~ 7rnob^DSc{0)p^. 

The value of 5*c(0) could then be obtained, at least in principle, comparing the amplitude of the interference pattern 
with the measured transmitted intensity. 



VI. SUMMARY AND OUTLOOK 



The example of neutron interaction with matter has been discussed inside the approach outlined in to describe 
the subdynamics of a microsystem interacting with a system having many degrees of freedom. The formal scheme 
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leads to a generator for the irreversible time evolution of the Lindblad form, whose expression relies on suitable 
choices for the potential term related to the T-matrix and the statistical operator describing the thermodynamic 
state of the system. In the example considered the main ingredient is given by the Fermi pseudopotential adopted 
to describe the neutron- nucleus interaction in impulse approximation. Then we obtain from (2.7), neglecting the 
incoherent contribution, the equation used by Sears to describe all neutron optical phenomena, as well as known 
expressions for the index of refraction. The incoherent contribution is necessary to fulfill the optical theorem and take 
diffuse scattering, that attenuates the coherent beam, into account. We have also shown how it may be connected to 
properties of the macrosystem, as expressed by the dynamic structure function. Furthermore possible experimental 
implications have been discussed in Sec. V. 

Even though it introduces a smaller correction the incoherent contribution is very important from the theoretical 
point of view. We expect that it will help studying the tricky borderline betwee n a p ure optical wavelike behavior 
and the fully incoherent particlelike one, based on a diffusion equation: in fact (2.10) leads in a direct way to the 
theory of Brownian motion , as stressed under (^|^); a somewhat similar treatment, in the case of an ideal gas, has 
been given in . It is not surprising that the incoherent contribute to the dynamics has grown out of a thoroughly 
quantum mechanical treatment, as shown by the typical quantum structure of the Lindblad equation, relying on 
non-commutating operators, in which an essential role is played by the statistical operator g, rather then by the wave 
function tj}. This point is of central relevance, since the terms which describe the incoherent dynamics cannot be 
introduced in the formalism of the wave function and are therefore unavoidably absent in an optical- like treatment, 
simply reminiscent of classical optical descriptions. 

We hope that this study of the emergence of incoherence in neutron-matter interaction will lead to a better under- 
standing of the general problem of irreversibility and of description of non-equilibrium systems. Typically coexistence 
of an incoherent particlelike behaviour, described by a Quantum Boltzmann equation, and a wavefunction description 
by means of Gross-Pitaevskii equation, is important for understanding Bose-Einstein condensation . In it is 
shown how the formalism we have used in the present paper copes with the more general problem of non-equilibrium 
macroscopic systems. However a systematic treatment of irreversibility in the very similar problem of atomic inter- 
ferometry involves QED and is a future challenge. 
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